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摘 要 : 将 生态 系统 服务 评估 利用 到 区 域 生态 保护 决策 中 的 重要 途径 之 一 是 生态 系统 服务 制图 。 本 研究 对 衡阳 


市 多 期 遥感 影像 进行 解 译 ， 
值 损 益 ; 并 在 ArcGIS 10.0 
的 时 空 变化 进行 制图 模拟 ， 


采用 当量 因子 法 估算 其 生态 系统 服务 价值 (ESV) 及 由 土地 利用 格局 转变 而 引起 的 价 
软件 支持 下 ,利用 空间 自 相 关 、 高 低 聚 类 及 重心 分 析 等 空间 统计 方法 ， 对 衡阳 市 ESV 
以 期 为 衡阳 市 的 生态 保护 规划 与 管理 提供 决策 辅助 。 研究 结果 表明 : 1)2000 一 2013 年 


T 


衡阳 市 的 ESV 略 有 降低 ， 


H 2000 年 的 2410 亿 元 降 至 2013 年 的 2394 亿 元 ; 转移 矩阵 分 析 发 现 林 地 及 耕地 转 为 


建设 用 地 ， 以 及 林地 转 为 耕 
失 区 主要 分 布 于 市 区 及 其 东部 区 域 ， 而 增值 区 主要 分 布 于 市 区 西部 区 域 .3) 衡 阳 市 ESV 


35; 2000—2013 年 ESV 损 


"n 


地 是 引起 衡阳 市 ESV 降低 的 最 主要 原因 。2) 衡 阳 市 的 ESV 高 值 区 对 低 值 区 呈 包 围 态 


的 空间 自 相关 与 高 低 值 聚集 现象 明显 ， 其 程度 趋 于 增强 ; ESV 重心 向 西 偏 北 9 方向 迁移 了 245.73 m. 可 考虑 通过 


限制 林地 及 耕地 向 建设 用 地 的 转变 、 实 施 “ 存 量 " 或 “ 减 量 " 用 地 政策 、 维 护 和 强化 区 域 整体 生态 系统 格局 的 连续 


性 等 措施 来 提高 研究 区 的 生态 系统 服务 价值 。 
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Abstract: Using ecosystem services assessment in regional ecological protection is one of the most important ways in ecosys- 


tem services mapping. In order to perfect the method of ecosystem services mapping, spatio-temporal variations in ecosystem 


services value (ESV) in Hengya 


ng Region was simulated based on interpretation of multi-period remote-sensing images of 


Hengyang City, estimation of ESV with the equivalent factor method and incurred profit-loss due to transformation of land use. 


Also the results of methods such 


as spatial autocorrelation, high-low clustering and barycenter analysis were analyzed in Ar- 


cGIS 10.0 environment. The results showed that: 1) the proportion of ESV components were different in the study area and 
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decreased with the sequence of biodiversity conservation, water conservation, soil formation and conservation, climate regula- 
tion, gas exchange, raw material, waste disposal, entertainment and food production. The ESV for Hengyang City decreased 
slightly from 241.0 billion Yuan to 239.4 billion Yuan during 2000-2013. The transformation of forests and cultivated lands 
into construction lands and the transformation of forests into cultivated lands were identified as the main reasons for the de- 
crease in ESV revealed by conversion matrix analysis. 2) Areas with low ESV were nearly encircled by those with high ESV. 
While regions with losses in ESV were mainly distributed around the urban and county districts, value-added regions were 
located in the western zone of the study area. 3) Spatial autocorrelation and value aggregation phenomena were obvious, which 
tends strengthened. The barycenter of ESV moved 245.73 meters west-northwest 9?. To improve ESV in the study area, meas- 
ures including limiting transformation of forests and cultivated lands into construction lands, conducting "stock" or “reduc- 
tion" policy on land use, and preserving and intensifying continuity of whole ecosystem patterns were required. Compared 
with previous studies, the spatial distribution of ESV was visualized in this study as more precisely expressing the income and 


transfer of ESV. This was critical for the promotion of more reasonable ecological protection planning and management works 


in the study area and beyond. 


Keywords: Ecosystem service value; Land use; Spatial statistical analysis; Hengyang City 


生态 系统 服务 是 指 自然 生态 系统 及 其 物种 所 


提供 的 能 够 满足 和 维持 人 类 生活 需要 的 条 件 和 过 
程 5”]。 对 其 价值 进行 定量 评估 有 利于 高 效 、 合 理 配 


置 竞 争 性 环境 资源 B， 有 利于 科学 制定 生态 环境 建 
设 补偿 政策 由。Costanza 等 最 先 对 生态 系统 服务 价 
值 (ecosystem service values, 以 下 均 用 “ESV”* 表 示 ) 
进行 了 定量 评估 ， 并 成 为 国际 可 持续 发 展 研究 的 热 
点 之 一 。 欧 阳 志 云 等 "1 侍 伯 杰 等 对 生态 系统 服 
务 价值 的 概念 、 效 用 及 评估 方法 进行 了 较 详 细 的 论 
xh; 谢 高 地 等 0 多 次 对 Constanza 的 研究 方法 进 
行 了 改进 , 得 出 了 适合 我 国 国情 的 ESV 当量 系数 ， 
利用 当量 因子 法 对 区 域 ES V. 进行 评估 已 被 国内 外 学 
者 广泛 采纳 。 此 外 ,综合 生态 参数 和 价值 当量 ,以 及 
从 生态 资产 角度 对 区 域 ESV 进行 评估 的 方法 也 得 到 
了 较 好 应 用 与 发 展 " M, 

土地 利用 /覆盖 变化 (LUCC) 影 响 生 态 系统 的 结 
构 和 功能 051， 其 对 生态 系统 维持 服务 功能 起 决定 
性 作用 59, 因此 ， 诸 多 学 者 利用 不 同方 法 在 不 同时 
间 、 空 间 尺 度 上 对 土地 利用 变化 引起 的 ESV 变化 
进行 了 研究 ("SI。 随 着 3S 技术 的 发 展 ,， 国内 外 已 
有 不 少 学 者 将 遥感 数据 与 生理 生态 模型 相 结合 ， 采 
用 当量 因子 进行 ESV 评估 ,并 对 其 空间 分 布 进行 
制图 模拟 ， 其 结果 为 区 域 保护 决策 与 管理 规划 的 制 
定 提供 了 很 好 的 支持 "” (BRIBI 3S 技术 进行 区 
域 ESV 损益 情况 的 空间 分 布 模拟 以 及 利用 空间 统 
计 分 析 来 揭示 ESV 的 时 空格 局 变化 规律 与 特征 的 
研究 尚 显 不 足 。 

本 研究 以 湖南 省 衡阳 市 为 例 ， 在 ArcGIS 10.0 软 
件 支持 下 ， 基 于 衡阳 市 多 年 遥感 影像 数据 及 经 济 发 
展 统计 数据 ， 对 其 ESV 进行 评估 ,进行 ESV 及 其 损 
益 的 空间 分 布 制图 模拟 ; 并 利用 空间 统计 分 析 方 法 ， 


探究 衡阳 市 ESV 的 时 空 动态 演化 特征 。 研 究 结果 可 
为 衡阳 市 进行 区 域 生态 安全 管理 及 发 展 生态 型 经 济 
提供 必要 参考 。 


1 研究 区 概况 

衡阳 市 位 于 湖南 省 中 南部 ， 湘 江 中 游 ， 衡山 之 
南 ， 是 湖南 省 域 副 中 心 城市 、 中 南 地 区 重要 工业 城 
市 、 区 域 性 物流 中 心 , 现 辖 5 KARK, EEK, 
HERK, RWE BEES 县 (衡阳 县 、 衡 南 县 、 衡 
LE, HWRE, MRE) RE 2? 个 县 级 市 (来 阳 市 、 
常宁 市 )。 地 处 东经 110932'16"—113?1632", 3b 2b 
26°07'05”~27°28'24”， 总 面积 约 15 310 km’ 近年 来 ， 
衡阳 市 经 济 与 社会 发 展 迅速 , 区域 土 地 利用 类 型 及 
其 空间 结构 也 随 之 发 生 了 较 大 变化 ， 进 而 影响 到 区 
域内 生态 系统 服务 价值 。 
2 研究 方法 
21 数据 来 源 与 处 理 

采用 衡阳 地 区 2000 年 、2005 年 、2010 年 及 
2013 年 同期 的 LandSatTM 遥感 影像 ， 空 间 分 辨 率 
30 mx30 mo 根据 研究 所 需 ,， 采 用 目 视 解 译 方式 ， 将 
研究 区 土地 利用 类 型 划分 为 耕地 、 林 地 、 草 地 、 湿地 、 
水 域 、 建 设 用 地 及 未 利用 地 共 7 个 类 型 得 到 研究 区 土 
地 利用 类 型 (图 1)， 其 总 体 精度 分 别 为 85.2%、86.5%、 
88.7%、87.4%，Kappa 值 超过 0.8， 达 到 本 研究 所 
需 。 相 关 经 济 与 社会 发 展 数据 来 源 于 衡阳 市 统计 
年 鉴 。 
2.2 ”研究 方法 
2.2.1 土地 利用 格局 转变 的 ESV 核算 

采用 基于 单位 面积 价值 当量 因子 的 方法 对 衡 
阳 市 土地 利用 格局 转变 引起 的 ESV 变化 进行 核算 。 
ESV 当量 因子 是 生态 系统 产生 服务 的 相对 贡献 率 。 
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1 衡阳 市 2000 £F. 2005 4E. 2010 年 和 2013 年 土地 利用 类 型 
Fig. 1 Maps of land use types in Hengyang City in 2000, 2005, 2010 and 2013 


根据 谢 高 地 等 2 的 定义 ， 中国 地 区 1 个 ESV 当量 因 
子 的 经 济 价值 等 于 当年 全 国平 均 粮 食 单产 市 场 价格 
的 1/7。 为 使 研究 结果 体现 区 域 特征 ,在 利用 研究 区 粮 
食 产量 进行 ESV 当量 因子 计算 的 同时 , 引入 基于 生物 
量 的 调整 因子 00 对 研究 区 的 ESV 价值 系数 进行 修正 。 
计算 研究 区 的 土地 利用 转移 矩阵 ， 并 核算 研究 
区 土地 利用 格局 转变 引起 的 ESV 损益 。 传 统 的 基于 
各 用 地 类 型 面积 的 ESV 分 析 模 型 ， 只 考虑 监测 期 间 
的 ESV 数值 的 变化 ,无 法 深层 次 地 反映 ESV 损益 的 
去 向 及 来 源 。 基 于 土地 转移 矩阵 分 析 ESV, 能 够 同 
时 考虑 各 用 地 类 型 的 ESV 的 空间 转移 分 布 和 数量 变 
化 方向 ,能 更 深入 揭示 ESV 的 变化 特征 。 
核算 过 程 所 涉及 到 的 计算 公 \ 式 见 表 1。 
2.2.2 ESV 演变 的 空间 统计 分 析 
地 理 数 据 受 空间 相互 作用 及 空间 扩散 的 影响 ， 
某 位 置 上 的 数据 与 其 他 位 置 上 的 数据 彼此 间 存 在 相 
互 依赖 的 现象 即 空 间 自 相关 性 。 对 ESV 在 一 定 地 
理 范围 内 不 同 空间 单元 间 的 相互 影响 和 依赖 程度 进 


行 分 析 , 是 把 握 其 演变 特征 与 规律 的 关键 。 采 用 Global 
Moran's I, Getis-Ord General G 探测 研究 区 ESV 在 特 
定 尺度 上 的 空间 自 相关 性 、 高 / 低 聚 类 。 Global 
Moran's I 用 于 分 析 ESV 的 整体 空间 自 相 关 性 ， 
Getis-Ord General G 用 于 探测 ESV 是 否 存 在 空间 集 
聚 。 重 心 分 析 是 一 种 探讨 地 理 现象 空间 分 异 及 其 动 
态 变 化 过 程 的 方法 ,运用 重心 概念 , 分 析 ESV 移动 


方向 及 距离 来 反映 区 域 ESV 的 空间 差异 , 探寻 其 动 
态 演变 过 程 。 

分 析 中 所 涉及 的 公式 见 表 2。 
3 结果 与 分 析 


3.1 ESV 的 变化 
3.1.1 ESV 结构 的 变化 

基于 研究 区 遥感 解 译 数 据 (图 1), 利用 上 述 ESV fh 
算 模 型 ， 得 出 衡阳 市 2000 年 、2005 年 、2010 年 及 2013 
{F ESV 值 ( 表 3), 并 计算 这 4 个 时 期 的 ESV 均值 , 制作 
衡阳 市 2000—2013 年 ESV 平均 分 布 格局 图 (图 2)。 
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表 1 生态 系统 服务 价值 核算 中 利用 到 的 公式 及 其 变量 
Table 1 Formulas and variables in the ecosystem services value (ESV) accounting 
公式 合 义 公式 内 容 变量 变量 合 义 Enn 
Meaning of formula Content of formula Variable Meaning of variable bere 
生态 系统 服务 价值 计算 公式 ESV = Yu xVO) ESV 生态 系统 服务 价值 Ecosystem services value Y-hm 7a! 
Formula of ESV n E 
Ak 土地 利用 类 型 的 面积 Area of land use type k hm* 
土地 利用 类 型 的 生态 系统 价值 系数 -2-1 
VC, Wi ¥:hm ^a 
Coefficient of ESV of land use type k 
生态 系统 第 /项 服务 功能 ESV; = $ (4ix VC) ESV 生态 系统 第 /项 服务 功能 价值 xx 
价值 计算 公式 f Value of ecosystem service function f 
Formula of ESV of item f VC; 土地 利用 类 型 上 的 第 /项 服务 功能 价值 系数 NEUE 
: Value of ecosystem service function f of land use type k 
生态 系统 服务 价值 系数 公式 _1 paxpyxB; Im VC 生态 系统 服务 价值 系数 Coefficient of ESV Y-hm ^a^ 
Formula of coefficient of ESV 7 f cu p 7A 
: 研究 未 年 全 国平 均 粮食 价格 2 
National average price of grain at the end of the period 8 
ō 研究 时 段 内 的 平均 粮食 年 产量 i 
Average annual grain yield during the study period £ 
中 国 陆地 第 ;类 生态 系统 的 第 7 项 服务 功能 的 当量 因子 5 
Di Value equivalent factor of ecosystem services function j in E 
Chinese terrestrial ecosystem i 
B. 第 i 类 生态 系统 的 生物 量 因 子 调整 系数 ("1 
Adjustment coefficient of biomass factor in ecosystem i 
土地 利用 格局 转变 的 ESV VC, = X(A, xVC, - VC,) eb 土地 利用 格局 转变 的 ESV 损益 系数 dpm 
损益 系数 计算 公式 ° ESV income coefficient of land use pattern transfer 
Coefficient formula of land VC, 土地 利用 类 型 上 转变 前 的 ESV 系数 Yhm?a 
use pattern transfer income ESV coefficient before transfer of land use type k 
表 2 生态 系统 服务 价值 空间 统计 分 析 中 利用 到 的 公式 及 其 变量 
Table 2 Formulas and variables in the spatial statistical analysis of ecosystem services value (ESV) 
ANTE SN 恋 量 首位 
a 公式 内 容 变量 变量 合 义 pan 
Ordinal Ha Content of formula Variable Meaning of variable Do 
formula variable 
I 空间 自 相 关 分 析 s pi in (x -x)(x;-x) I 全 局 空间 自 相关 指数 Index of spatial autocorrelation 
计算 公式 -A Mo B l . 
Formula of spatial Aa” PEU Xi, Xj 空间 单元 六 /的 ESV 指数 ESV index of space unit ij — Y-hm 7a! 
autocorrelation x ESV 的 均值 Average value of ESV Yhm^a' 
W; ZHAN EEE Binary space weight matrix — 
n 空间 单元 个 数 Number of units 一 
2 高 / 低 聚 类 指数 计 — o eR S [vv G — G 指 数 取 值 范围 为 [0, 1]， 值 越 高 越 趋向 于 高 值 聚 集 。 
算 公 式 Ma C fija CEC), B. Z(G) 显 著 时 , 检测 出 高 值 聚集 ; G 值 <E(G)， 
Formula of (为 ) 且 Z(G) 显 著 时 , 检测 出 低 值 聚 集 (E(G) 为 G 的 理论 期 望 
high/low cluster f&, Z(G)7J G 的 理论 方差 值 。The G index range was [0, 
index ; i ; 
1], the higher value tends to higher value aggregation. = 
When G > E(G) and Z(G) significantly, the high value 
aggregation detected. When G < E(G) and Z(G) signifi- 
cantly, the low value aggregation detected. E(G) is 
theoretical expectation of G, Z(G) is theoretical value of 
variance of G. 
3 冷 / 热 点 指数 计算 G;* GEREI Z 得 分 越 高 , 热点 聚集 越 紧密 ,表明 单元 周围 的 属 
公式 
Formula of [ Jum E x-lw | 
hot/cold index i n Ls 
The higher the Z score of the G;* index is, the closer the hot spot E 
aggregation, which shows that the attribute value around unit is 
relatively high. And: X I x S= 
n ia 
Xj 空间 单元 7 的 ESV 指数 变化 量 z 
Variation of ESV index of spatial unit j 
E 重心 计算 公式 > J xX, y 表示 “重心 ”的 坐标 Coordinates of the barycenter 一 
Formula of bary- yu K F-a " 
center on’ n Xi, yi 要 素 i 的 坐标 Coordinates of the factor i = 
n 要 素 总 数 Number of factor 一 
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ANTEA SN 75 E yr 
AE oe 公式 内 容 变量 变量 合 义 E 
Ordinal seen DAE Content of formula Variable Meaning of variable e 
mula variable 
5 重心 移动 距离 计 D, -dq, -Yy +(X, -Xy] 1/2 Dn m 年 份 到 n 年份 间 重 心 移动 的 距离 im 
Box Moving distance of barycenter from year m to year n 
Formula of bary- QC, ym) (Xm yo) AQ, y) 表 示 第 mm 年 和 第 年 的 区 域 重心 坐标 
center moving (Xn Yn) (Xm, Ym), (Xn, Yn) represent regional barycentric coordi- m 
nates in year m and n. 
C 常数 , 取 111.111， 是 把 地 全 标 单位 (经 纬度 ) 换 算 成 平 
面 距离 (km) 的 系数 
Constant, value is 111.111, the coefficient converting 
the whole standard unit (latitude and longitude) into the 
plane distance (km). 
表 3 衡阳 市 2000—2013 年 生态 系统 服务 价值 结构 变化 
Table 3 Structure changes of ecosystem services value of Hengyang City in 2000— 2013 
服务 价值 Service value («105 Y) 2000— 2005 2005— 2010 2010— 2013 2000— 2013 
入 功能 
É 服务 功能 . 价值 变化 ”变化 率 ”价值 变化 变化 率 价值 变化 变化 率 ”价值 变化 ”变化 率 
t 
gor cq 2000 2005 2010 2013 Value Change Value Change Value Change Value Change 
change rate change rate change rate change rate 
CD — (9) COW (Q0) Ð W — (103) — (9) 
气体 交换 
331.29 330.76 | 330.17 329.20 -0.53 —0.16 —0.60 —0.18 —0.96 —0.29 —2.09 —0.63 
Gas regulation 
气候 调节 
. 气候 调节 : 330.89 331.49 330.49 329.58 0.61 0.18 —1.00 —0.30 —0.91 —0.27 —1.30 —0.39 
Climate regulation 
水 源 涵养 . 360.94 361.83 359.68 358.79 0.89 0.25 -2.15 —0.59 —0.89 —0.25 -2.15 —0.60 
Water conservation 
土壤 形成 与 保护 
Soil formation and 346.78 346.13 344.85 343.85 -0.65 —0.19 一 1.28 —0.37 —1.00 —0.29 —2.93 —0.85 
protection 
^ 
废物 处 理 219.58 220.77 218.37 217.89 1.18 0.54 -2.40 —1.09 —0.48 —0.22 —1.70 —0.78 
Waste treatment 
生物 多 样 性 保护 
Biodiversity 365.61 365.18 364.08 363.04 —0.43 —0.12 一 1.11 —0.30 一 1.04 —0.29 -2.57 —0.70 
conservation 
食物 生产 71.18 71.01 70.19 67.00 -0.17 —0.24 —0.82 -1.15 —0.19 —0.28 -1.18 -1.66 
Food production 
原材料 
. 222.43 221.94 221.74 221.09 -0.49 —0.22 —0.20 —0.09 —0.65 —0.29 —1.34 —0.60 
Raw material 
娱乐 休闲 
Entertainment 161.40 161.61 161.07 160.63 0.20 0.12 —0.54 —0.34 —0.43 —0.27 —0.78 —0.48 
culture 
合计 Total 2410.10 2410.72 2400.632394.07 . 0.65 0.03 —10.09 —0.42 —6.56 —0.27 —16.04 —0.67 


计算 结果 表明 : OMA ESV 价值 构成 方面 来 看 , 衡 
阳 市 各 单项 生态 功能 价值 所 占 比 例 的 大 小 依次 为 生物 
多 样 性 保护 > 水 源 洱 养 > 土壤 形成 与 保护 > 气候 调节 > 
气体 交换 > 原材料 > 废物 处 理 > 休 闲 娱 乐 > 食 物 生 产 。 前 
5 项 服务 功能 的 价值 量 相对 较 高 且 相 差 不 大 ，2000 一 
2013 年 平均 占 总 功能 价值 量 的 比例 分 别 为 15.16%、 
14.99%、14.37%、13.75%、13.74%; 原材料 与 废物 处 
理 的 功能 价值 量 次 之 , 这 2 项 各 占 总 功能 价值 量 的 
9.23%、9.12%; 食物 生产 与 娱乐 休闲 的 价值 量 最 低 , 两 
者 各 占 总 功能 价值 量 的 2.94% 和 6.71%。@) 从 ESV 变化 
情况 来 看 , 2000 一 2013 年 衡阳 市 ESV 变化 较 缓 , 所 有 单 
项 ESV 均 略 有 降低 , 总 量 减 少 0.67%。 在 2000 一 2005 
年 间 , 衡阳 市 ESV 量 增加 了 0.62 亿 元 , 是 由 于 其 生态 服 
务 功 能 中 的 废物 处 理 、 水 源 洱 养 、 气 候 调 节 及 娱乐 休闲 
的 价值 量 略 有 增加 , 分 别 增 加 了 0.54%、0.25%、0.18% 


和 0.12%; 2005—2010 年 间 及 2010—2013 年 间 , 所 有 
单项 ESV 均 呈 减少 趋势 。 

ESV 平均 分 布 格局 图 显示 : 衡阳 市 ESV 高 值 区 
大 体 分 布 于 靠近 市 界线 的 山林 地 ， 而 低 值 区 多 分 布 
于 市 域 范围 的 较 中 心 处 ,此 范围 内 地 势 较 缓 ， 适 宜 
城镇 建设 及 农耕 开发 。 整 体 看 来 ， 高 值 区 与 低 值 区 
均 具 有 较 好 的 整体 性 及 连续 性 ， 高 值 区 对 低 值 区 呈 
包围 态势 , 低 值 区 以 市 区 为 核 向 各 县 城区 心 蔓延 ， 
与 衡阳 市 的 盆地 地 形 暗合 。 
3.1.2 土地 利用 格局 转变 的 ESV 损益 

根据 图 1 计算 衡阳 市 土地 利用 格局 的 转移 和 矩阵， 
并 结合 ESV 计算 公式 ， 得 出 由 于 土地 利用 格局 转变 
引起 的 ESV 损益 结果 ( 表 4)， 并 利用 ArcGIS 软件 中 
的 栅 格 计算 器 及 插值 分 析 得 出 各 研究 时 段 的 ESV 损 

空间 分 布 图 (图 3)。 
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ESV (X) 
Bl o s54.s31 5 
B 72 854831 51~110 129.396 4 
BE 110 129 396 5-154 181.155 

[7] 154 181.155 1-198 232.913 6 
[| 198 232.913 7-432 046.093 8 


0 15 30 60 —— 高 速 公路 Expressway 
km 一 一 一 铁路 Railway 


2 2000—2013 年 衡阳 市 生态 系统 服务 价值 (EVS) 平 均 空间 分 布 格局 
Fig.2 Spatial distribution pattern of average ecosystem services value (EVS) in Hengyang City from 2000 to 2013 

计算 结果 表明 : (02000—2005 4, BFE ”时段 更 为 剧烈 。 增 值 中 贡献 较 多 的 土地 利用 转变 
地 转变 引起 ESV 增值 5.13 亿 元 , 损失 4.51 亿 元 ， ”方式 为 耕地 向 林地 、 水 域 的 转变 ,高达 70908918 
共计 增值 0.62 亿 元 。 增 值 中 贡献 较 多 的 土地 利用 ” 值 来 自 于 这 两 种 转变 方式 ; 而 林地 转 为 耕地 及 建 
转变 方式 为 耕地 转 林 地 、 耕 地 转 水 域 、 林 地 转 水 ，” 设 用 地 是 引起 此 时 段 ESV 损失 的 主要 原因 ， 总 损 
域 以 及 水 域 转 湿地 这 4 种 , 约 86% 的 增值 来 自 于 失 量 的 87% 由 这 两 种 转变 引起 。(3)2010 一 2013 年 
这 4 种 转变 方式 ; 耕地 、 林 地 转 为 建设 用 地 是 引起 i, ESV 增值 5.85 亿 元 , 损失 12.41 亿 元 ,共计 损 
此 时 段 ESV 损失 的 主要 原因 ,总 损失 量 的 73% 由 K 6.57 亿 元 。 其 中 耕地 与 林地 的 相互 转变 引起 的 
这 两 种 转变 引起 。@2005 一 2010 年 间 , ESV 增值 ESV 损益 基本 持平 ， 造 成 ESV 损失 的 主要 原因 在 
56.27 亿 元 , 损失 66.37 亿 元 ， 共 计 损 失 10.09 亿 元 ， ”于 林地 及 耕地 向 建设 用 地 的 转变 ,这 两 种 转变 引 
显示 出 此 时 段 内 的 土地 格局 变化 较 另 外 两 个 研究 E ESV 损失 达 6.1 亿 元 。 


人 


2000 一 2005 年 gi 


ll 


e 


n 


增益 Profit 0 1530 60 km "^ -铁路 Railway 
- 损失 Loss 高 速 公路 Expressway 
3 2000—2013 年 衡阳 市 生态 系统 服务 价值 (ESV) 损 益 分 布 空间 格 局 
Fig.3 Spatial distribution patterns of profit and loss of ecosystem services value (ESV) in Hengyang City from 2000 to 2013 
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ESV 损益 空间 分 布 格局 图 显示 : (002000—2005 
年 间 , ESV 损失 区 主要 分 布 于 市 区 、 衡 东 县 城区 北部 
以 及 来 阳 市 城区 南 稍 偏 东部 处 ; 增值 区 主要 分 布 于 
衡阳 县 东部 、 邢 东 县 西部 、 衡 南 县 西南 部 以 及 常宁 
市 东部 区 域 。@2005 一 2010 年 间 ，ESY 的 损益 分 布 
较 另 外 两 个 研究 时 段 复 杂 , 损失 区 主要 分 布 于 衡阳 
县 东部 与 市 区 相 邻 区 域 、 衡 南 县 城区 西南 部 区 域 、 
来 阳 市 南端 的 铁路 及 高 速 公路 周边 区 域 、 衡 东山 县 
的 南部 及 北部 区 域 、 衡 山 县 南部 以 及 南岳 区 西 稍 偏 
北部 区 域 ; 增值 区 主要 分 布 于 衡 南 县 东北 部 与 市 区 
相 接 处 、 衡 阳 县 西北 部 、 姥 东 县 西北 部 、 衡 东 县 东 
部 、 常 宁 市 西北 部 以 及 来 阳 市 东北 及 东南 部 区 域 。@) 
2010—2013 年 间 , ESV 损失 区 主要 分 布 于 衡 南 县 东 
部 、 衡 东 县 南部 以 及 常宁 市 东南 部 区 域 , 另外 市 区 以 
及 衡阳 县 、 衡 东 县 、 衡 南 县 与 市 区 接壤 区 域 的 ESV 
也 损失 较 明 显 ; 增值 区 主要 分 布 在 衡阳 县 北部 、 衡 
山 县 西北 部 、 衡 阳 县 西部 、 姥 东 县 南部 以 及 来 阳 市 
南部 区 域 。 


32 ESV 的 空间 自 相关 与 高 / 低 聚 类 分 析 

计算 结果 ( 表 5) 显 示 : 空间 自 相 关 分 析 中 , 4 个 
年 份 的 7 值 较 高 , 同时 引入 标准 化 统计 量 Z(7D) 用 以 检验 
空间 自 相 关 的 显著 性 水 平 ，2Z(7)= [7-E(7)]/ 
VVar(7) [其 中 ECD、Var(D) 分 别 为 指数 7 的 理论 期 户 
值 与 理论 方差 值 ]。 结 果 显 示 , Z(D) 值 全 部 高 于 1% 显 
著 性 水 平 标准 值 ， 这 表明 衡阳 市 ESV 的 空间 分 布 具 
有 和 较 强 的 空间 自 相 关 性 ， 空 间 分 布 格局 呈 聚 集 模 
式 ; 另外 4 个 时 期 的 7 指数 与 Z(D) 值 呈 先 下 降 后 上 
升 的 趋势 ， 表 明 衡 阳 市 的 ESV 空间 自 相 关 性 先 趋 
向 减弱 再 趋 于 增强 。@ 高 / 低 聚 类 分 析 中 , 4 个 年 份 
ESV 的 G 指数 无 明显 大 于 E(O) 值 , 但 Z(G) 值 较 显 
Z, 表明 衡阳 市 ESV 的 空间 分 布 存在 高 值 及 低 值 
聚集 模式 ; Z(G) 值 呈 先 下 降 后 上 升 趋势 ， 表明 ESV 
的 高 值 以 及 低 值 的 聚集 程度 在 2000 一 2005 年 间 有 
所 弱化 ， 其 后 又 一 直 增 强 并 超过 2000 年 时 的 聚集 
程度 。 


RS 2000—2013 年 衡阳 市 生态 系统 服务 价值 ESV) 空 间 自 相关 值 (D) 与 高 / 低 聚 类 指数 
Table 5 Spatial autocorrelation & high/low clustering index of ecosystem services value (ESV) in Hengyang City in 2000-2013 


ESV 空间 自 相 关 Spatial autocorrelation of ESV 


ESV 的 G 统计 Getis-Ord general G of ESV 


指数 index 指数 index 
2000 2005 2010 2013 2000 2005 2010 2013 
Moran's 7 0.157 1 0.156 0 0.166 0 0.166 6 G ]E-6 ]E-6 ]E-6 ]E-6 
Z() 8.997 0 8.933 5 9.503 2 9.537 6 E(G) ]E-6 ]E-6 ]E-6 ]E-6 
\ == Z(G) 5.400 0 5.321 2 5.676 1 5.674 4 


3.3 ESV 高 值 区 重心 迁移 轨迹 

根据 表 2 中 的 公式 4， 计 算 4 个 年 份 的 ESV 加 
RED, 制作 ESV 重心 迁移 轨迹 图 (图 4)， 并 统计 其 
迁移 距离 (图 5)。 由 图 4 可 看 出 , 2000— 2013 年 衡阳 市 
的 ESV 重心 一 直 落 入 于 珠 晖 区 长 尾 冲 附近 。 图 5 显示 
2000—2013 年 衡阳 市 的 ESV 重心 迁移 非常 缓慢 ,主要 
向 西 稍 偏 北方 向 移动 , 大体 上 呈 “U”* 型 态势 。 ESV 重心 
迁移 具有 一 定 的 阶段 性 , 经历 了 “缓慢 一 剧烈 一 缓慢 ” 
的 变化 阶段 ，2000 一 2013 年 间 的 总 迁移 直线 距离 约 为 
247 m, 分 阶段 来 看 , 2005—2010 年间 ESV 重心 平均 年 
际 移动 距离 最 长 , 为 ESV 重心 迁移 的 剧烈 阶段 ， 其 余 
两 个 阶段 平均 移动 距离 相差 不 大 。 


4 讨论 与 结论 

本 研究 以 当量 因子 法 为 基础 ， 对 衡阳 市 的 遥感 
影像 进行 解 译 , 得 到 7 种 类 型 的 土地 利用 方式 ， 利 
用 生物 量 指标 对 ESV 当量 因子 进行 修正 ,估算 各 年 
份 研 究 区 的 ESV, 并 对 研究 区 的 ESV 分 布 及 损益 的 
空间 分 布 进行 制图 模拟 ， 运用 空间 自 相 关 、 高 低 聚 


类 及 加 权重 心 求 算 等 空间 统计 分 析 方 法 ,探究 衡阳 
T ESV 的 时 空 动态 演化 特征 。 研 究 结果 可 为 衡阳 市 
决策 者 制定 符合 客观 实际 的 土地 利用 结构 、 进 行 环 
境 保 护 与 管理 生态 系统 提供 参考 , 促进 衡阳 市 经 
济 、 社 会 与 生态 的 协调 发 展 。 

研究 发 现 : 002000— 2005 年 , 衡阳 市 ESV 总 量 
略 有 上 升 ， 这 应 该 主要 得 益 于 湖南 省 自 2000 年 开展 
实施 的 退耕 还 林 工 程 ; 而 2005 一 2013 年 衡阳 市 ESV 
总 量 的 持续 下 降 ， 可 能 是 衡阳 市 城市 化 进程 中 侵占 
林地 及 耕地 带 来 的 必然 问题 。@ 衡 阳 市 ESV 的 平均 
分 布 格局 为 外 围 高 内 心 低 ， 这 主要 与 衡阳 市 的 盆地 
地 形 有 关 , 已 有 不 少 研究 发 现 地 形 因 子 对 区 域 的 生 
态 环境 质量 影响 显著 1。 影响 衡阳 市 ESV 损益 
变化 的 主要 原因 为 林地 、 耕 地 及 建设 用 地 三 者 间 的 
交互 转变 , 损益 分 析 图 显示 ,，ESV 损失 区 主要 发 生 
在 市 区 及 各 县 城区 所 在 地 ， 这 应 是 城市 建设 使 大 量 
耕地 及 林地 转 为 建设 用 地 ， 卫星 影 像 图 显示 ,衡阳 
市 的 建设 用 地 面积 2000—2013 年 共 增 加 了 
18 504.72 pm。 衡阳 市 的 ESV 重心 在 向 西 稍 偏 北 
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4 2000—2013 年 衡阳 市 生态 系统 服务 价值 (ESV) 重 心 落 入 处 ( 左 图 ) 及 其 局 部 放大 图 ( 右 图 ) 
Fig.4 Location of ecosystem services value (ESV) gravity centers in Hengyang (left picture) and its partial enlarged detail (right 
picture) from 2000 to 2013 
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5 2000—2013 年 衡阳 市 生态 系统 服务 价值 (ESV) 重 心 迁 移 轨迹 
Fig. 5 Moving track of ecosystem services value (ESV) gravity centers in Hengyang City from 2000 to 2013 


方向 移动 。 这 可 能 与 其 各 区 、 县 的 经 济 发 展 与 人 口 
增长 的 不 平衡 有 关 。 衡 阳 市 统计 年 鉴 “ 显示 , 衡阳 
市 的 经 济 及 人 口 存在 空间 分 异 , 来 阳 市 、 常 宁 市 及 
衡 南 县 的 经 济 发 展 与 人 口 增长 规模 较 祁 东 县 、 衡 东 
县 及 衡山 县 速度 快 , 衡阳 市 的 县 域 经 济 发 展 水 平 趋 
势 主 要 为 东 高 西 低 ， 这 和 与 本 研究 的 ESV 重心 变化 轨 
迹 相 呼 应 。 

需要 指出 的 是 : ONEI ESV 估算 方法 得 到 的 
结果 很 可 能 存在 较 大 差异 , 目前 ESV 评估 方法 可 大 
致 分 为 “基于 单位 服务 功能 价格 的 方法 "和 “基于 单 
位 面积 价值 当量 因子 的 方法 ”两 类 ,相对 于 服务 价格 
法 ， 当 量 因子 法 由 于 直观 且 所 需 数 据 量 少 ,， 且 该 方 
法 特别 适用 于 区 域 及 全 球 尺度 的 ESV 评估 ， 便 于 进 
ITAA GIS 进行 分 析 ， 因 此 研究 采用 此 方法 。 凶 当 


量 因子 法 存在 缺陷 之 处 ， 如 研究 区 域内 的 某 种 用 地 
类 型 的 面积 未 发 生变 化 , 但 其 景观 指数 如 斑 块 数目 、 
平均 板块 面积 、 聚集 度 等 发 生 了 变化 , 这 势必 会 对 ESV 
产生 影响 , 虽然 生物 量 法 对 此 能 够 有 所 修正 , 但 不 能 完 
美 解决 此 问题 , 因此 有 学 者 指出 应 进行 不 同 空间 层次 、 

不 同 尺度 的 ESV 评估 实证 分 析 "…"”。@@) 除 用 生物 量 法 
进行 修正 外 ,有 学 者 还 根据 研究 区 农田 生态 系统 的 实 
际 情况 , 对 其 食物 生产 的 ESV 系数 进行 修正 59， 或 利 
用 突变 理论 对 我 国 的 陆地 生态 系统 的 单位 面积 ESV 系 
数 进行 修正 "1 不 同 的 区 域 修正 方法 必 将 产生 不 同 的 
ESV 最 终 评估 结果 。(4) 若 能 细 分 土地 利用 类 型 及 分 析 
更 高 分 辨 率 栅 格 , 研究 结果 将 更 为 可 信 ; 另外 限于 篇 幅 , 
本 文 主 要 着 重 于 ESV 的 平均 分 布 、 损 益 进行 制图 ， 以 
及 空间 统计 方面 的 相关 分 析 , 未 对 ESV 的 时 空 分 异 特 
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征 及 其 变化 的 驱动 机 制 进行 分 析 , 这 将 在 以 后 的 研究 
中 给 予 完 善 。 
本 研究 主要 结论 如 下 : 


1)2000 一 2013 年 , 衡阳 市 的 ESV 由 2000 年 的 
2 410 亿 元 降 至 2013 年 的 2 394 亿 元 ; 研究 区 中 生物 
多 样 性 保护 、 水 源 涵养 、 土 壤 形 成 与 保护 、 气 候 调 
节 及 气体 交换 这 5 项 服务 功能 较为 突出 , 食物 生产 与 
娱乐 休闲 的 价值 量 最 低 。 衡阳 市 的 ESV 高 值 区 对 低 值 
区 呈 包 围 态势 , 高 值 区 主要 分 布 于 市 界 四 周 的 山林 地 , 
中 心 区 域 由 于 地 势 较 缓 利于 农耕 及 建设 ,因而 成 为 
ESV 的 低 值 区 , 这 与 衡阳 市 的 盆地 地 形 暗 合 

2) 林 地 及 耕地 转 为 建设 用 地 ， 以 及 林地 转 为 耕 
地 是 引起 衡阳 市 ESV 降低 的 最 主要 原因 ; ESV 损失 
区 主要 分 布 于 市 区 、 衡 南 县 、 衡 阳 县 东部 区 域 、 来 
阳 市 城区 南部 以 及 常宁 市 城区 东南 部 。 

3) 衡 阳 市 的 ESV 表现 出 明显 的 空间 自 相 关 与 高 
低 值 聚集 现象 , 2000 一 2013 年 总 ESV 呈 降 低 趋势 ， 


但 其 自 相关 与 聚集 程度 趋 于 增强 ， 即 其 空间 分 布 的 
破碎 化 程度 在 降低 ; ESV 重心 一 直 落 入 于 市 区 南部 


稍 偏 西 与 衡 南 县 车 江 镇 的 交接 处 ， 且 迁移 速度 非常 
缓慢 ， 主 要 向 西 稍 偏 北 方向 移动 ， 呈 “缓慢 一 剧烈 一 
缓慢 ”的 趋势 。 

根据 研究 结论 ， 可 考虑 以 下 措施 实现 对 衡阳 市 
生态 系统 服务 的 有 效 管控 : 

1) 限 制 林地 及 耕地 向 建设 用 地 的 转变 ; 鼓励 植 
树 造 林 、 退 耕 还 湖 ， 以 抵消 建设 用 地 扩张 而 带 来 的 
ESYV 损失 。 

2) 对 ESV 损失 区 重点 监测 ， 实 施 “ 存 量 ” 或 “ 减 量 ” 
用 地 政策 以 严格 控制 其 范围 ,引导 ESV 系数 低 的 土 
地 利用 类 型 向 高 的 类 型 转变 。 

3) 在 编制 区 域 土地 利用 总 体 规划 中 ， 可 通过 维 
护 和 强化 区 域 整体 生态 系统 格局 的 连续 性 ， 加 强 对 
林地 、 湿 地 及 水 域 的 保护 ,以 提高 研究 区 ESV。 
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